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Abstract 

A nonlinear model of the quantum harmonic oscillator on two-dimensional spaces 
of constant curvature is exactly solved. This model depends of a parameter A that 
is related with the curvature of the space. Firstly the relation with other approaches 
is discussed and then the classical system is quantized by analyzing the symmetries 
of the metric (Killing vectors), obtaining a A-dependent invariant measure dfj-x and 
expressing the Hamiltonian as a function of the Noether momenta. In the second part 
the quantum superintegrability of the Hamiltonian and the multiple separability of the 
Schrodinger equation is studied. Two A-dependent Sturm-Liouville problems, related 
with two different A-deformations of the Hermite equation, are obtained. This leads 
to the study of two A-dependent families of orthogonal polynomials both related with 
the Hermite polynomials. Finally the wave functions ^m.n and the energies Em^n of 
the bound states are exactly obtained in both the sphere S'^ and the hyperbolic plane 
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1 Introduction 



This article can be considered as a sequel or continuation of a previous paper [1] which was devoted 
to the study of a quantum exactly solvable one-dimensional nonlinear oscillator with quasi-harmonic 
behaviour. Now, our idea is to extend the results and present a similar analysis but for the quantum 
version of the two-dimensional nonlinear system. We follow the approach of [1], which contains the 
fundamental ideas and motivation, and we will also make use of some properties discussed, at the 
classical level, in [2] (other related papers are [3]-[6]). 

The following nonlinear differential equation 

(1 -I- \x^) X — \xx^ + X = , (A a constant) , (1) 

was study in Refs. [7, 8] as an example of a nonlinear oscillator. In Lagrangian terms, the equation 
(1) can be obtained from the following function 




that clearly displays two very interesting characteristics: the potential ^(A) has a nonpolynomial 
character (this is not a problem of an harmonic oscillator perturbed by higher order terms of the 
form Ax^"* with m > 1), and the kinetic term depends of a position-dependent mass. 

Let us briefly comment these two characteristics. 

The form of the potential V{\) is shown in Figures I and II for several values of A (A < in 

Figure I and A > in Figure II). We see that for A < the potential is a well with a boundless wall 
at = 1/|A| and for A > we have that V ^ (1/2)(q;^/A) for x— > ± oo. It can be proved that 

1. If A < then the general solution of (1) is given by 

X = Asui{uj t + (f)) , = {1 + \A^)u? . 

2. If A > then the general solution is given by 

X = Asm.{uj t + <p) , = (1 + A^^)a;^ , 
when the energy E is smaller than the value E^^x = a^/(2A), and by 

x = Bsinh(m-F0i), = (AS^ - l)a;2, 
when the energy E is greater than E^ x- 

The Schrodingcr equation involving the potential A(,x^/(1 + gx^)) has been studied by many 
authors making use of different approaches [9]- [30]. In some cases the idea was to study the 
Hamiltonian H = — (f /dx^+x^+X{x^ / (l+gx^)) by applying perturbative, variational or numerical 
techniques previously used for the system = — cP/dx^ +x'^ + Ax^"* with m > 1. It is important 
to note that in most of these papers the derivative part of the Schrodinger equation was the standard 
one, that is, the equation arising from a classical Hamiltonian with a quadratic term of the form 
(l/2)p^ and leading to a derivative term of the form — d^/dx^, or to the corresponding two or 
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three-dimensional versions involving the Laplace operator in or E^. Nevertheless, we point out 
that Mathews and Lakshmanan studied in Ref. [10] the following quantum Hamiltonian 



H = \[\{p\{l-9x')}^ + 



(1-5x2) J' 



where the notation {A, 5}+ = AB + BA is used. 

The second important point is the presence of a position-dependent mass m = {1 + \x^) ^ since, 
if the mass becomes a spatial function, then the quantum version of the mass no longer commutes 
with the momentum. Therefore, different forms of presenting the kinetic term in the Hamiltonian H, 
that are equivalent at the classical level, lead to different and nonequivalent Schrodinger equations 
[31]- [42]. This is an old question that remains as an important open problem in the theory of 
quantization. 

In spite of these two characteristics, the quantum version of this A-dependent non-linear oscillator 
has been proved to be exactly solvable [1] . The question of the order ambiguity in the quantization 
of the Hamiltonian was solved by introducing a prescription obtained from the analysis of the 
properties of the classical system (existence of a Killing vector and a A-dependent invariant measure) 
and concerning the problems arising from the nonpolynomial character of the potential, all of them 
disappear when ^(A) is studied with the appropriate quantization of the position-dependent kinetic 
term. 

On the other hand the following two-dimensional Lagrangian 



was proposed in Ref. [2] as the appropriate two-dimensional generalization, at the classical level, 
of the Lagrangian (2) . In fact the general solution of the Euler-Lagrange equations which are given 




) 



9 2\/ \2l^/^ \ 2 



2 I 2 

X +y , 



by 



{I + Xr"^) X - \ [x^ + ip' + \{xy 
{l + \r^)y-\[x^ + f + \{xy 



yx)"^ ] .X + .T = , 
yxf ]y + y = Q , 



is: 



1. If A < then the general solution is given by 



X = Asm.(u)t + 4)i) , y 



B sm.(u)t + (p2) , 



for all the values of the energy E. 



2. If A > then the general solution is given by 



X = Asm{ijjt + 4>i) , y 



B sin(a;t -|- ^2) , 



when the energy E is smaller than a certain value E^^^x, and by 



X = A sinh(r2i + 4)i) ., y = B sinh(Oi -|- ^2) 



when the energy E is greater than this value, E > Eo,,\- 
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In both cases the coefficients A and B are related with the coefficient a and the frequency lo 
(oscillatory motions) or with a and $7 (unbounded motions). So we have "quasi-harmonic" nonlinear 

oscillations in the case of bounded motions and high energy scattering solutions when A > 0. 
Moreover, the analysis of this nonlinear system proved the existence of a relation with the linear 
harmonic oscillator on the sphere or on the hyperbolic plane with the parameter A playing 
the role of the (negative of the) curvature k. 

The main objective of this article is to quantize this two-dimensional nonlinear oscillator as a 
deformation of the harmonic oscillator in the sense that 

1. All the fundamental properties of the linear system continue to hold for A 7^ but modified 
in a A-dependent way. 

2. The limit when A ^ is well defined and when A = all the characteristics of the quantum 
harmonic oscillator are recovered. 

The idea is to prove that it is exactly solvable and to obtain the energies and the corresponding 
wave functions. Some of the main questions to be discussed in this paper can be summarized in 
the following four points: 

• Relation of this A-dependent nonlinear model with the harmonic oscillator on spaces of con- 
stant curvature. 

This quasi-harmonic nonlinear oscillator is related, at the classical level, with the harmonic 
oscillator on the three spaces of constant curvature {S^, IR^, H^). Now this relation is 
considered for the quantized systems. 

• Analysis of the transition from the classical A-dependent system to the quantum one. 

The two-dimensional A-dependent kinetic term possesses three Noether symmetries. The 
main idea is to quantize the system by using as Hilbert space the space L'^{'\R,dfi\) where 
diJ,x is a measure invariant under the Killing vectors associated to the Noether symmetries. 

• Schrodinger separability and 'quantum superintegrability'. 

The two-dimensional Schrodinger equation is not separable in {x,y) coordinates because of 

the A-dependent coupling between the two degrees of freedom; nevertheless it is proved that it 
admits separability in several coordinate systems. The existence of this multiple separability 
is a property related with 'quantum superintegrability'. 

• Exact resolution of the A-dependent Schrodinger equation and families of A-dependent or- 
thogonal polynomials. 

Two A-dependent Sturm-Liouville problems related with two different A-deformations of the 
Hermite equation are obtained. This leads to the study of two A-dependent families of 
orthogonal polynomials both related with the Hermite polynomials. 

In more detail, the plan of the article is as follows: In Sec. 2 we study the relation of this 
nonlinear model with the harmonic oscillator on spaces of constant curvature. In Sec. 3 we first 
introduce a A-dependent measure, we use it for introducing a quantization rule and we obtain the 
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A-dependent Schrodinger equation. Sec. 4 is devoted to the analysis of the A-dependent Schrodinger 
separabihty and to solve two Hermite-related equations and Sec. 5 to obtain the eigenfunctions 

^m,n and energies Em^n- In Sec. 6 we briefly analyze the existence of another possible approach 
and its relation with the presence of the angular momentum. Finally, in Sec. 7 we make some final 
comments. 



2 On the relation of this nonhnear model with the har- 
monic oscillator on spaces of constant curvature 



Although the first studies of the harmonic oscillator on curved spaces are rather old (the last chapter 
of Ref. [43], that is called 'Nichteuklidische Mechanik', is devoted to the dynamics on spaces with 
curvature; it first studies general properties and then it consider the harmonic oscillator as a 
particular case; the approach is mainly Newtonian), Ref. [44] is usually considered as the more 
relevant paper for the modern approach to this system (by modern we mean that it studies subjects 
such as dynamical symmetries or quantum dynamics). We recall that the harmonic oscillator is a 
system that is well defined not only in the Euclidean plane but also in the other two-dimensional 
spaces of constant curvature, sphere S"^ and hyperbolic plane H^. 

In differential geometric terms, the three spaces with constant curvature, sphere S"^ {k > 0), 
Euclidean plane E^, and hyperbolic plane (k < 0), can be considered as three different situations 
inside a family of Riemannian manifolds = (S'^,E^,i?^) with the curvature k as a parameter 
At G IR. In order to obtain mathematical expressions valid for all the values of k, it is convenient 
to make use of the following /c-trigonometric functions 



cos if At > 0, 

Ck(x) = -j 1 if K = 0, 

, cosh-v/— K X if K < 0, 



S«(x) = < 



sm^/Kx 



■■ sinh-v/^ 



KX 



if At > 0, 
if K = 0, 
if K < 0, 



then the expression of the differential element of distance in geodesic polar coordinates {R, $) on 
the family = (S^, E^, iJ^), can be written as follows 



dsl = dR^ + Sl{R) d^^, 



so it reduces to 



dsf = dR^ + (sin^ R) d^^ , ds^ = dR^" + R^ d^^ , dst^ = dR^" + (sinh^ R) d^^ , 

in the three particular cases of the unit sphere, the Euclidean plane, and the 'unit' Lobachewski 
plane. Note that R denotes the distance along a geodesic on the manifold M^; for example, in the 
spherical k > case, R is the distance of the point to the origin (e.g., the North pole) along a 
maximum circle. 

If we make use of this formalism then the Lagrangian of the harmonic oscillator on is given 
by [45, 46] 

IL(«) = ih {vr + Sl{R)vl) - (i) a' tUR). (3) 
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where the K-dependent tangent is defined in the natural way Tk{R) = Sk{R)/ Ck{R)- In this way, 
the potential of the harmonic oscillator on the unit sphere, on the Euclidean plane, or on the unit 
Lobachewski plane, arise as the following three particular cases 

Ui{R) = {^) a'^ tan^ R , UoiR) = i^) R^ , C/_i(i?) = (^) tanh^ . 

The Euclidean oscillator Uo{R) appears in this formalism as a parabolic curve making a separation 
between two different situations (see Figure III). Note also that in spherical k > case, this 
Lagrangian describes in fact two oscillators with centers in the north {R = 0) and south {y/nR = tt) 
poles and with a boundary barrier in the equatorial circle. 

Next we study the behaviour of 1L(«;) under two different changes of variables. 



1. Let us consider the K-dependent change (i?, $) (r, (j)) given by 

r = SK{R), = ^, \ = -K. 
Then the Lagrangian ]L(/«) becomes 



2 V l + Ar 



, 2 2 



a 

T 



Therefore, if we change to Cartesian coordinates (x, y) we arrive to 

^^^^ = \ (tTaT^) + + ^ (^^^^ - y'-A - T (tTaT^) ' 



2 2,2 



This function is just the Lagrangian obtained in Ref. [2] as the natural generalization of 
the one-dimensional Lagrangian L{x,Vx'-iX) for the nonlinear equation (1) of Mathews and 
Lakshmanan. 

2. Let us now consider the K-dependent change {R, $) (r', (p) given by 

r' = T«(i2), = 
Then the Lagrangian ]L(«;) becomes 



Lh{k) 



1 



.,'2 



^ 2 12 

-ar . 
2 



2 \^(1 + Kr'2)2 (l + Kr'2) ^ 
Therefore, if we change to Cartesian coordinates {x,y) we arrive to 



1 



1 



2 (1 + Kr'2)2 



^ 2 12 

-a r . 
2 



„'2 



X +y , 



This function is the Lagrangian studied by Higgs in Ref. [44] (the study of Higgs was originally 
limited to a spherical geometry but the idea can be easily extended to the hyperbolic space). 
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We note that these two changes are correct and both radial variables, r and r', are well defined. 
In the hyperbolic k < case the two functions Sk{R) and Tk{R) are positive for i? > and 
concerning the spherical k > case this property is also true because then R is restricted to a 
bounded interval. 

The situation can be summarized as follows. We have obtained three alternative ways of describ- 
ing the harmonic oscillator on spaces of constant curvature: the original system IL(«;) and the two 
other approaches, L(A) and Lh{k), obtained from it. Of course, everyone of these three different 
approaches, IL(k), L(A) and Lh{k), has its own characteristics and advantages. 

A two-dimensional manifold M can be described by using diff'erent coordinate systems. If we 
consider it as an imbedded submanifold of IR^, then the points of M can be characterized by 
the three external coordinates, as e.g. {x,y,z), plus an additional constraint. Nevertheless, in 
differential geometric terms, a more appropriate approach is to develop the study by using two- 
dimensional systems of coordinates intrinsically defined in M (and without make reference to the 
external space). The Lagrangian TL{k) is directly defined on the manifold = (S*^, E^, H"^) and it 
uses the expression of the differential element of distance ds^ in geodesic polar coordinates (R, 
(see the Appendix). Therefore, in differential geometric terms, this approach can be considered as 
more formally correct than the other two. 

The Higgs approach [44, 47] consider the motion on S", embedded in the Euclidean space E""*"^, 
by means of a central (also known as gnomonic) projection on a plane fl" tangent to at a 
chosen point. This particular formalism leads to a dynamics that is described, when n = 2, by the 
Lagrangian Lh{k). This approach is very interesting because it states a direct relation between 
the motion on a curved space, the sphere S"", and the motion on a plane. In fact, it has been 
studied by many authors (see e.g. [48]-[59] and references therein) mainly in relation of the theory 
of dynamical symmetries. 

The A-dependent Lagrangian L{\) has a certain similarity with the Lagrangian of Higgs but 
nevertheless it does not coincide with it: in the model of Higgs k (or A) is present in the kinetic 
term T in a different way and the potential V appears as ^-independent; wc will see that this 
affects, via the Legendre transformation, to the Hamiltonian formalism. As wc have seen in Sec. 
I (Introduction) one of the advantages of this approach is that the Euler-Lagrange equations can 
be directly solved and the general solution has a rather simple form that can be interpreted as 
"quasi-harmonic" nonlinear oscillations. It is clear that this A-dependent formalism seems very 
appropriate for solving equations or for other related calculus. 

In what follows we will focuss our attention on the Hamiltonian dynamics determined by the 
A-dependent Lagrangian L{\). 



3 Quantization and A-dependent Schrodinger equation 



Let us start our study considering the following Lagrangian 



= \ (t 



+ Ar2 



a 
~2 



1 + Ar2 



r'' = x'' + y\ (4) 



where the parameter A can take both positive and negative values; of course it is clear that for 
A < 0, A = — |A|, the function (and the associated dynamics) will have a singularity at 1— |A| = 
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and we shall restrict the study of the dynamics to the interior of the interval < 1/|A| where the 
kinetic energy function is positive definite. 

The Legendre transformation is given by 

Vx — XJy Vy + XJx 



Px 



1 + Ar2 ' 1 + Ar2 ' 

and the expression of the A-dependent Hamiltonian turns out to be 



1 



= 7J Px+Py + >' [xPx + yPyY 



+ 



a 



2 Vi + Ar2 



(5) 



The transition from the classical system to the quantum one is a difficult problem because of the 
ambiguities in the order of positions and momenta. In Refs. [1, 3] the one-dimensional nonlinear 
oscillator was quantized by using a prescription obtained from the existence of a one-dimensional 
Killing vector X{\) and a A-dependent measure dii\ in M preserved by -'^(A). Next we prove that 
this approach, that was successful for the one-dimensional system, admits a direct generalization 
to this more difficult two-dimensional case. We must begin with an analysis of the symmetries of 
the kinetic energy term. 

It was shown in [2] that function T(A) representing the kinetic energy 



is invariant under the action of the vector fields Xi(A), X2(A), and Xj, given by 



Xi(A) 

X2{\) 

Xj 



vT+Ar2 — , 
ox 

Vl + Ar2 

dy 
d__ d_ 

dy ^ dx ^ 



in the sense that, if we denote by Xf., r = 1,2, J, the natural lift to the tangent bundle (phase 
space IR^xlR^) of the vector field Xr, 

d 



XliX) 



Vl + Ar2 / dvx 
d 



/ i , X 2 ^ , \ (X}!x±y3i\ o 



dy 



dy 

d_ ^ 

dx ^ dv,. 



_d_ 

dVcr 



then the Lie derivatives of T(A) with respect to Xl{\) vanish, that is 

X*(A)(r(A)) =0, r = l,2,J. 

It is known that a symmetric bilinear form in the velocities (v^, Vy) can be considered as associated 
to a two-dimensional metric ds^ in IR^. In this particular case, the function r(A) considered as a 
bilinear form determines the following A-dependent metric 



ds'^{\) 



Vl + Ar2 



(H- A y^) dx^ (1-1- A x^) dy"^ - 2A xy dx dy 



(6) 
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Thus, in differential geometric terms, the three vector fields Xi{X), X2{X), and Xj, must be con- 
sidered as three Killing vector fields (infinitesimal generators of isometries) of ds^(A). These three 
symmetries of the kinetic term determine three associate Noether momenta (a Noether momentum 
is a constant of motion for the geodesic motion) given by 

that become 

-Pi (A) = a/ 1 + Ar^ , P2{\) = + Xr'^ Py, J = xpy-ypx, 

in the Hamiltonian formalism. We note that the form of the angular momentum J is preserved by 
the Legendre map, in the sense that we have xpy — ypx = xvy — yvx (this is another one of the 
differences with the Higgs model). 

Proposition 1 The only measure on the space that is invariant under the action of the three 
vector fields Xi(X), X2{X), and Xj, is given by 



Vl + Ar- 
up to a constant factor. 

Proof: The most general expression for a volume two-form on the space IR^ is given by 

Lo = p{x, y) dx A dy 
Then the Lie derivatives of uj under Xi(A), X2(A), and Xj, are given by 

Lxt,oj = (Vl + Xr^ dx A dy + p {dVl + Xr^) A dy 
Lx-^u) = (^Vl + Xr'^ dx A dy + dx A p {dVl + X r^) 

The condition Cxjd^x = implies that p{x,y) must be a function /(r) of r. Then the two other 
conditions, Cxidjjix = and jCx2dP'X = 0, lead to (for r 7^ 0): 

/— 1 df Xf 
V1 + Xr^ - — + =0, 
r dr V 1 + A 

with general solution given by 

Vl + Ar2 

where k is an arbitrary numerical constant. g 

We will consider this proposition as the fundamental point for the study of transition from the 
classical system to the quantum one. In fact, this property suggests us to work with functions and 
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linear operators defined on the space obtained by considering the two-dimensional real plane IR^ 
endowed with the measure d^x given by 



(7) 



This means, in the first place, that the operators Px an Py representing the quantum version of 
of the Noether momenta momenta Pi an P2 must be self-adjoint not in the standard space L^(IR) 
but in the space {JR, d/ix) . If we assume 

d 



Py = -in Vl + Ar2 



dx ' 
d_ 

dy ■ 



then we arrive to the following correspondence 



as well as 



d_ 

^ dy ^ dx 



d 



dy 



^ dy ^ ) ' 



in such a way that the quantum version of the classical Hamiltonian 

+~ 

that can be rewritten as follows 



(8) 



IS 



H = 



2rn 

+ A — 
2m 



(1 + Ar^)^ + A.- 



2m 



52 



92 



92 



5x dy 



(l + Ar^ 
A_ _ 



52 ^ 



d_ 
dy 



(9) 



The first important property of this Hamiltonian is that it admits the following decomposition 

H = Hi+H2- XJ^ 
where the three partial operators Hi, H2 y are given by 



Hi = 
H2 = 



2m 
2m 



(1 + Ar 
(l + Ar 



9x2 
2^ 5^ 



+ Ax 



;i 2 + Ay 



d_ 

dx 
d_ 

dy 



+ 



1 + Ar2 

y2 



+ Ar^ 
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in such a way that the the total Hamiltonian H commutes, for any value of the parameter A, with 
each one of the three partial terms 

[H,Hi] = 0, [H,H2]=0, [H,J']=0. 

The vanishing of these three commutators means that the A-dcpcndent Hamiltonian (9) describes 
a quantum superintegrable system [60]- [71]. This property was analyzed at the classical level in 
Ref . [2] ; now we see that the quantization rule we have applied preserves the superintegrability. 

Now, if we consider the Schrodinger equation 

as we have the following property 

[i?i,i?2-Aj2] =0, [Hi - X.'P ,H2]=0, [Hi + H2,J^]=0, 

then, we have three different sets of compatible observables and therefore three different ways of 
obtaining a Hilbert basis of common eigenstates. 

1. The two operators Hi and H2 — \J^ are a (complete) set of commuting observables; therefore 
they represent two quantities that can be simultaneously measured. Thus, the first way of 
looking for ^' is as a solution of the following two equations 

^1 * = El * , (^2 - A * = E2j * . 

In this case the total energy is given hy E = Ei+ E2j and the associated wave function can 
be denoted by •^{Ei,E2j). 

2. The two operators Hi — XJ^ and H2 are a (complete) set of commuting observables. Thus, 
the second way of looking for is as a solution of the following two equations 

(Fi - A P) * = Eij * , ^2 ^' = ^2 ^ • 

In this case we have E = Eij + E2 and ^ can be denoted by ^{Eij,E2). 

3. The third (complete) set of commuting observables is provided by Hi + H2 and J^. So in 
this case we have 

{Hi + F2) ^ = £^12 * , J ^ = j ^ . 

Thus, the two physically measurable quantities are E12 and the angular momentum j, the 
total energy is given by £^ = £^12 — Aj^ and the wave function so defined can be denoted by 

^{El2,j). 

The existence of these three alternative descriptions arises from the presence of the term A J 
inside the kinetic part of the Hamiltonian. Notice that the second approach is just the symmetric 
of the first one but, although they are closely related, they lead however to different solutions with 
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different properties; that is, "^{Ei,E2j) / ^'(E'lj, £"2). This fact is a consequence of the nonhnear 
character of the model since in the hnear hmit, when X ^ 0, then both descriptions coincide. 



Let us consider the following quantum Hamiltonian 



H = - 

+y 



2m 

,2 



52 



^2 



(1 + ^^')^ + ^^ 



dx"^ ^ dx dy 



dx J 2m 
d__ d_ 

dx ^ dy 



dy^ ^ dy 



+ A 



2m 



dy^ 



+ (l)^a(„ + AA)(- 



(10) 



where we have slightly modified the value of the parameter g that now is given hy g = ma^ + Xha 

(this is done to coincide with the notation of the one-dimensional system in [1]). It is also convenient 
to simplify this function H by introducing adimensional variables {x, y, A, e) defined by 



( —) 
V v ma ) 



y 



A 



(— )a, E = in, 



a) e . 



f2 



in such a way that the following relation holds 

l + Ar^ = 1 + Af2 
The Schrodinger equation takes then the following form 

92 . , a 1 . 1 



~2 I -2 



2 52 , ^2 d 

+ y TTTiT - 2xy 



_ _ 

2 

. d 



(1 + A^=)^ + A!/^ 



92 



5x dy 



dx 



y 



dy 



f2 



(11) 



4 Resolution of the Schrodinger equation I 



In the following all the variables, parameters and equations are adimensional. Nevertheless, and 
for ease of notation, we will drop the use of the tilde and write the variables just as x, y, r and so 
on. 



4 . 1 Separability 

The Schrodinger equation (11) is not separable in Cartesian {x,y) coordinates because of the A- 
dependent term. Nevertheless, at this point we recall that it was proved in [2] that, at the classical 
level, the Hamilton-Jacobi equation 

fc) + +t( i + a(/+,2) )=2£- 

admits separability in the following three different orthogonal coordinate systems: 
1. A-dependent coordinates {zx,y) with Zx defined by Zx = x j \J \ -\- h.y^ . 
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2. A-dependent coordinates {x,Zy) with Zy defined by Zy = y/vT+~Ax^. 

3. Polar coordinates {r,(p). 

The expression of the potential V{X) in these three systems is 

= -(^—) 

1 / 

~ 2 VI + Ar^ 



1 + Az2 

+ 



1 + Az2 



Since it is known the existence of a close relation between (additive) classical Hamilton-Jacobi 
separability and (multiplicative) quantum Schi^odinger separability [72, 73] it seems natural to make 
use of these three coordinate systems for the study of this quantum problem. 

Next we start our study with the first coordinate system. 

Using {zx,y) coordinates, the Schrodinger equation becomes 



'l + Az^l^ 



/L-tn^\ u-_ / Az^ X u 
Vl + AysJ dz^ \l + Ay^J dz^ 



1 

2 



& 



d 



(12) 



We assume a solution in the form 



^{z,y) = Z{z,)Y{y), 



where Z and Y are, respectively, functions of Zx and y alone. Substituting this expression in the 
Eq. (12) we arrive to 



1 1 

2 Z 



1 1 

2 Y 



(1 + A zl) Z" + (A Zx) Z']+{\){l + k) {yT)^) 
:i + A y-'f Y" + 2A y (1 + A y2) y ] _ (1) (i + A) + (i + Ay^) g 



that leads to the two following ordinary equations 



1 1 

2 Z 

1 1 

2 Y 



{l + Kzl)Z" + {Kzx)Z']+{\){l + K) [y^^zI 

(l + Ay2)2y" + 2Ay(l + Ay2)r']+(^)(i + A)y2_(i + ^y2)g^_^^ 



where [i denotes the separation constant. 

Consequently the A-dependent Schrodinger equation is in fact separable in the [zx-, y) coordinates. 
Thus the two-dimensional problem has been decoupled in two one-dimensional equations. 
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4.2 Power series resolution of the Z-equation 

The first equation to be solved is 

{l + kzl)Z" + {kz,)Z'-{l + k) {y^^) Z = -2^Z (13) 



This equation coincides (up to the appropriate changes of notation) with the equation of the 
one-dimensional nonlinear oscillator studied in Ref. [1]. Therefore, we directly explain the charac- 
teristics of the solution. 

Firstly, using the following factorization for the function Z 

Z(^,,A)=p(z,,A)(l + Az2)-i/(2A), (14) 

the function -p = p{zx, A) must satisfy the differential equation 

(1 + A zl)p" + (A - 2)zxp' + (2/x -l)p = 0, (15) 

that represents a A-deformation of the Hermite equation. Secondly, this new equation can be solved 
by the use of power series expansions. Assuming 

oo 

p{zx, A) = J2 = co(A) + ci(A) Zx + C2(A) zl + ... 

n=0 

the following A-dependent recursion relation is obtained 

Cn+2 = (-1) . [ n (A rz - 2) + (2/i - 1) 1 , n = 0, 1, 2, . . . 

[n + 2)[n + 1) L J 

The general solution p{zx), defined in the interval z^ < 1/A, is given by the linear combination 
p = C{jpi + cip2 where pi{zx) and P2izx) are the two solutions determined by pi{0) = 1, p'liO) = 
and P2(0) = 0, ^2(0) = 1, respectively. If there exists a certain integer m such that the coefficient 
IJ- = IJ-mis given by 

2;Um = 2 m + 1 — A , 
then we have Cm 7^ 0, Cm+2 = 0, and one of the two solutions is a polynomial of order m. 
The polynomial solutions are given by 



Even index (even power polynomials) 

r=p 

= E 



r=p 

2r 
C2rZx 



r=0 



«0 

[2-Ap'][2-A(p' + 2)][2-A(p' + 4)]...[2-A(p' + 2(r-l))], 



C2r = {-lY^p'{p'-2){p'-A)...{p'-2{r-l)) 



where we have introduced the notation p' = 2p. More specifically, the expressions of the first 
solution pi{zx), in the particular cases of p' = 0, 2, 4, are given by: 

Vo = 1, 

V2 = 1-2(1-A)z2, 

= l-4(l-2A)z2^(^)(l-2A)(l-3A)4-4A)4 
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• Odd index (odd power polynomials) 

r=p 
r=Q 

■ [2-A(p' + l)][2-A(p' + 3)]...[2-A(p' + (2r-l))], 

where we have introduced the notation p' = 2p + 1. More specifically, the expressions of the 
second solution P2{zx) for p' = 1, 3, 5, are given by: 

^1 = z^, 

n = z,-(^)(l-3A)4 + (^)(l-3A)(l-4A)4. 

4.3 Power series resolution of the F-equation 

The second equation to be solved is 

(l + A/)y" + (2Ay)y'- (l + A)(^^^)y + 2ey = 2(^^3^)y (16) 

that, although it has certain similarity with the Eq. (13), it docs not coincide with it (two differ- 
ences: the coefficient 2A and the //-dependent right-hand term). The main reason for this asym- 
metry is that, when introducing separability in the Schrodinger equation, the angular momentum 
term was displaced into this second equation. 

We start our study with the following two steps. 

Step 1. Introduction of a new quantum number 

It is convenient to decompose the energy e as the following sum 

e = iJ, + v 

where u is a new parameter. Then the equation (16) transforms into 

2 

(1 + Ay2)y" + (2Ay) y' - (1 + A - 2A/x) (yTA^) ^ + = 

which looks more similar to previous first equation for the function Z. 
Step 2. Factorization of the function Y 

Let us rewrite the previous equation as follows 

2 

(l + Ay2)y" + (2Ay)y'-G^(^-^^)y + 2z.y = 0, G2 =l + (l-2/.)A. 
Firstly, it can be verified that the function *oo defined by 

*oo = (l + AyV''''/^'^^ 
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satisfies the following property 



Thus, ^'oo is the exact solution of the Eq. (16) in the very particular case of u = (1/2)G^ and 
can be considered as representing, in the general case 21/ j/^ G^, the asymptotic behaviour of the 
solution. Consequently, this property suggests the following factorization 



y(y,A) = g(y,A)(l + Ay2)-G./{2A)^ 

and then the new function q{y, A) must satisfy the differential 

(1 + A + 2(A - + {2iy-G^)q = 0, 

that turns out to be a new A-deformation of the Hermite equation. 
Assuming a power expansion for the solution 

oo 

q{y, A) = ^ Cn{A) = co(A) + ci(A) y + C2(A) + . . . 

n=0 

the equation leads to 

oo 

^ [(n + 2) (n + 1) Cn+2 + A n(n - l)cn + 2(A - G^) n c„ + (2i/ - G^) c„] y'' 

n=0 

and we obtain the following A-dependent recursion relation 



(17) 
(18) 



C„+2 = (-1) 



(n + 2)(n + l) 



An (n - 1) - G^(2n + l) + 2u 



Note that this relation shows that, as in the A = case, even power coefficients are related among 
themselves and the same is true for odd power coefficients. In both cases, having in mind that 



lim, 





= hm„^oo 







An(n + 1) - G^(2n + 1) + 2z/ 
(n + 2)(n + l) 



the radius of convergence R is given by 



R 



Hence, when we consider the limit A — >^ 0, we recover the radius i? = oo of the Hermite's equation. 

The general solution is given by the linear combination q = cgqi + cig2 where qi{y) and q2{y) 
are the solutions determined by qi{0) = 1, q'i{0) = and q2{0) = 0, q2{0) = 1, respectively. In the 
very particular case of the coefficient v be given hy v = Vn with 

2i/n = (2n + 1)G^ — n(n + 1) A , (n is an integer number), 

then we have Cn^O, c„+2 = 0, and one of the two solutions becomes a polynomial of order n. 
The polynomial solutions are given by 
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Even index (even power polynomials) 

C2r = i~lY^p'ip'-2){p'-4)...{p'-2ir-l)) 

[2G- A{p' + 1) ][2G^ - Aip' + 3)][2G^- A{p' + 5) ] . . . [2G^ - A(y + 2r - 1) ] , 

where we have introduced the notation p' = 2p. More specifically, the expressions of the first 
solution qi{y), in the particular cases of p' = 0,2,4,6, are given by: 



Vo = I, 

Va = l-2(2G^-5A)y2 + (i)(2G^_5A)(2G^-7A)/. 



V2 = 1 - (2G^ - 3Ay , 



• Odd index (odd power polynomials) 

r=p 

'P2p+1 = "2r+l 2/^''"^^ 

r=0 

= (-l)'+i^-^(p'-l)(p'-3)...(p'-(2r-l)) 

■ [2G^ - A ip' + 2) ][2G^ - A{p' + 4) ] . . . [2G^ - A (p' + 2r) ] , 

where we have introduced the notation p' = 2p+l. More specifically, the expressions of the 
second solution q2{y) for p' = 1,3,5, are given by: 

Vi = y, 

V3 = y-(^)(2G^-5A)y3, 

V5 = y-(^)(2G^-7A)y3 + (^)(2G^-7A)(2G^-9A)/. 



5 Eigenfunctions ^rn,n and energies Em,n 

5.1 Sturm-Liouville problems and A-dependent Hermite polyno- 
mials 

We have obtained two different A-dependent deformations of the Hermite equation as well as the 
general solution and the particular polynomial solutions. Nevertheless in quantum mechanics the 
important point is not the equation by itself but the associate Sturm-Liouville problem. In this 

case we have two problems, one for the Z-equation and other for the "F-equation. Moreover each 
one of them splits into two: one for A < (spherical case) and other for A > (hyperbolic case). 

Spherical A < case The first A-dependent differential equation 

(1 + Azl)p" + (A - 2) za^p' + (2/i -l)p = (19) 
can be reduced to self-adjoint form by making use of the following integrating factor 

M..) = (l + A.^)-(A+2)/(2A)^ 
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in such a way that we arrive to the following expression 

d 
dzx 



A{zx,A)-^ +(2/x-l)r(^^,A)p = 0, (20) 



dzx 

where the two functions A = A{zx,A) and r = r{zx,A) are given by 

VTTA^ 1 



A 



(1 + Az2)l/A' V'l + A^^(l + Az2)l/A ' 

that, together with appropriate boundary conditions, constitute a Sturm-Liouville problem. 
If A is negative the problem is defined in the bounded interval [— a\ , a\] with ca = 1/ -v/IaJ- 



The function A{zx,A) vanishes in the two end points zi = — and Z2 = and the problem is 
singular. The boundary conditions prescribe that the solutions must be bounded functions at the 
two end points of the interval. The eigenvalues are the quantized values of the parameter /x, i.e. 
/j^m, m = 0, 1, 2, . . ., and the eigenfunctions the associated polynomial solutions. 

If A is positive the problem is singular since is defined in the whole real line IR. The solutions must 
decrease when ^ ± oo in such a way that their norms, determined with respect to the weight 
function r{zx), be finite. Therefore the eigenfunctions are again the A-dependent polynomials Vm, 
m = 0,1,2,... 

Proposition 2 The eigenfunctions Vm{zx-,A), m = 0,1,2, ... of the de Sturm-Liouville problem 
of the Eq. (20) are orthogonal with respect to the function r = (1 + A 2:2)~(^/2+iM) . 

Proof: This statement is just a consequence of the properties of the Sturm-Liouville problems. 
Because of this the polynomial solutions Vm-, m = 0, 1, 2, . . ., of the equation (20), satisfy 

/■"A V^{Zx,A)Vn{Zx,A) 

I ^ — „ — , dZx = 0, m n, A<0, 



p Vra{ZxA)Vn{ZxA) , ^ , . ^ . 

/ TT-m — ., dZx = 0, m n, A>0. 

7-00 l + KzlY/^Jl + Kzi ^ ' 



and 

In the A > the integral is defined on a infinite interval, the following property must be 

Scttisfi-Gci 

, hm zx WmizxA)? (1 + A^^)-(VA+i/2) =0. 

The consequence is that if A > then the quantum number m is limited by the condition m < 1/A, 
and there is only Ma eigenvalues and eigenfunctions where Ma denotes the greatest integer lower 
than 1/A. I 

The following "Rodrigues formula" 

Hmizx^A) = (-1)" W^/^+'/' ^ [ WT W- (i/^+V2) ] , w, = \^kzl, (21) 

leads to a family of A-dependent Hermite polynomials which are proportional to Vm 

'>im = kmVmizx^K) , m = 0, 1, 2, . . . 
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where km are constants that in the first cases are given by 

ko = fci = (2-A), 

= - (2 - 3A) , k3 = - 3(2 - 3A)(2 - 5A) , 

k4 = 3(2-5A)(2-7A), Aig = 15(2 - 5A)(2 - 7A)(2 - 9A) . 

Alternatively we can obtain these polynomials by using the following function 

Tiz^, A) = (l + A {2tz^ - t^)) ^'^ (22) 

as a generating function 

(l + A {2tz, -e))' =J2 (— r) '^miz,, A) , (23) 
where we have used the notation Hm for the coefficients of the Taylor series. We obtain 

'Hm = gmVm{Zx,J^) , m = 0, 1, 2, . . . 

where are constants that in the first cases are given by 

50 = 1 , 51 = 2 , 

92 = -2, 53 = -12(1 -A), 

54 = 12(1 -A), 55 = 120(1-A)(1-2A). 

We can define the A-dependent Hermite functions Zm by 

Z„(z^,A) =7^^(z^,A)(l + Az2)-i/(2A)^ m = 0,l,2,... 

then the above statement admits the following alternative form: The A-dependent Hermite func- 
tions Zjn{zx,A) axe orthogonal with respect to the weight function r = + Az'^: 

/OA ^ rllA dz.. 

Zm{zx,A)Zn{zx,A)r{zx,A)dzx= Zmizx, A) Zniz^, A)— =0, m^n, A<0, 

-OA J-aA y/l + Az^ 

and 

/oo ^ PCX) ^2 

Zm{zx,A) Zn{zx,A)r{zx,A)dzx = Zm{zx,A)Zn{zx,A)—=^== = 0, mj^n, A>0. 
-oo J-oo \Jl-\-Az^ 

Figures IV and V show the form of the function ^2(2^x5 ^) for several values of A (A < in Figure 
IV and A > in Figure V). 

Summarizing, the final solution of the Sturm-Liouville problem for the function Z(^Zx) is: 
• Spherical A < case: 

Zm{Zx,A) = nm{Zx,A) (1 - |A| ^2) V(2|A|) , 

fJ-m = {rn -\- -) -\- -rri^ \A\ , m = 0, 1, 2, . . . , m, . . . 
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• Hyperbolic A > case: 



A) = Urniz., A) (1 + A V(2A) ^ 



IJ-m = {m + -) - -m^ K , m = 0, 1,2, . . . ,Ma. 



The second A-dependent differential equation 



ao q" + aiq + a2q = Q , Gm = 1 — rnk , 

ao = l + Ay^ ai = 2(A-G^)y, a2 = 2v-G, 



(24) 



is not self-adjoint since u'q / ai but it can be reduced to self-adjoint form by making use of the 
following integrating factor 



that, together with appropriate conditions for the behaviour of the solutions at the end points, 
constitute a Sturm-Liouville problem. It is to be pointed out that the boundary conditions are 
in fact different according to the sign of A; therefore we arrive to, no just one, but two different 
Sturm-Liouville problems: 

If A is negative the range of the variable y is limited by the restriction < 1/|A|. In this case 
the problem, defined in the bounded interval [— ca, oa] with oa = l/-\/|A[, is singular because the 
function q{y, A) vanishes in the two end points yi = — a\ and y2 = a\. The conditions to be 
imposed in this case lead to prescribe that the solutions q{y, A) of the problem must be bounded 
functions at the two end points, yi = — a\ and y2 = a\, of the interval. It is clear that this leads 
to the above mentioned polynomial solutions. 

If A is positive the variable y is defined in the whole real line IR and, therefore, the Sturm- 
Liouville problem is singular. The solutions q{y, A) must be well defined in all JR, and the boundary 
conditions prescribe that the behaviour of these functions when y ^ ztoo must be such that their 
norms, determined with respect to the weight function r{y), be finite. It is clear that in this case 
the solutions of the problem are again the A-dependent polynomials Hm, m = 0, 1, 2, . . . 

Proposition 3 The eigenfunctions 7^„(y,m, A), n = 0, 1, 2, . . . of the de Sturm-Liouville problem 
of the Eq. (25) are orthogonal with respect to the function r = {1 + Aj/^)"*^""/^, = 1 — mA. 

Proof: This statement is just a consequence of the properties of the Sturm-Liouville problems. 
Because of this the polynomial solutions Vn = 'Pn{y,'m,^), n = 0,1,2,..., of the equation (25), 
satisfy 



ao 

in such a way that we arrive to the following expression 



^(y) = {—) e/(«i/»o)<i2/ = (1 + ^y2yGm/A ^ 



— [B{y, m. A) — J + {2u - Gm) r{y, m. A) g = , 
where the two functions B = B{y, m, A) and r = r(y, m, A) are given by 

B = (1 + A y2)l-G^/A ^ r = {l + A y2)-G^/A . 



(25) 




r s, A < 0, 
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and 



r Vriy,m,A)V.{y,m,A) 



(l+Ay2)G.n/A 



In the A > CcLSG ^ clS the integral is defined on a infinite interval, the following property must be 

satisfied 

, hm y [Vn{y, m, A)]^ (1 + A y2)m-l/A ^ ^ 

and as a consequence, the quantum number n is limited by the condition 

1 1 

n < - — m— - . 
A 2 

That is, for every value of m there is only A'a eigenvalues and eigenfunctions where N\ denotes 
the greatest integer number lower than 1/A — m— 1/2. g 

The "Rodrigues formula" for these new family of A-dependent Hermite polynomials is given by 
7^„(y,m,A) = (-irW«-/^^[w-;T^-«'"/A]^ Wy = l + Ky\ (26) 

in such a way that we obtain 

'Hn = knVn{y,rn,k), n = 0,1,2,... 

with the following values for the first constants 

feo = 1, fci = 2(G„, - A) , 

k2 = -2(G™-2A), fc3 = -12(G^-2A)(G„-3A), 

^4 = 12(G„^-2A)(G„-3A), A^s = 120(G^ - 3A)(G^ - 4A)(G„ - 5A) . 

The (A, m)-dependent function defined by 

Hy, t, A) = (l + A {2ty - e)) ^^-/^-'/'^ (27) 
is a generating function with the following power expansion 

(l + K{2ty-t')y ' ^(-)W„(y,m,A)r (28) 

m=0 

where we have used the notation Tin for the coefficients of the Taylor series. The first (A, m)- 
dependent Hermite polynomials obtained in such a way have the following expressions 

'Hn = 9nVn{y,m,A) , n = 0,1,2,... 

where the constants gi, i = 0,1,2, ... ,5, take the values 

go = I, gi = i2Gm - A) , 

52 = -(2G^-A), 53 = -3(2G„- A)(2G„-3A), 

54 = 3(2G„-A)(2G„-3A), 55 = 15(2G^ - A)(2G^ - 3A)(2G„ - 5A) . 
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The A-dependent Hermite functions ¥„, are defined by 

y„(y,m,A) =H„(y,m,A)(l + Ay2)-GW(2A)^ n = 0,1,2,... 

and the above statement admits the following alternative form: The A-dependent Hermite functions 
Yn{y,m, A) are orthogonal with respect to the weight function r = 1: 

/OA _ rOA 

Yr{y,m,A)Ys{y,m,A)r{y,A)dy = Yr{y,m,A)Ys{y,m,A)dy = 0, r s , A<0, 
- OA J — 0,f^ 

and 

/OO I'OO 
Yr{y,m,A)Ys{y,m,A)?{y,A)dy = Yr{y,m, A)Ys{y,m, A) dy = , r^s, A>0. 
-OO J —OO 

Summarizing, the final solution of the Sturm-Liouville problem for the function Y(y) is: 

• Spherical A < case: 

Yn{y, m, A) = m. A) (1 - | A| y^) Gr^/i^W) , G„, = 1 + m |A| , 

i'n = {n+ -)Gm + -^n{n + l)\A\, n = 0, 1, 2, . . . , m, . . . 

• Hyperbolic A > case: 

Yn{y, m, A) = H„(y, m. A) (1 + Ay^y ^-/(^A) ^ Gm = l-mA, 
i^n = {n+^)Gm-^n{n + l)A, n = 0,1,2, ... ,Na. 

5.2 Final solution 

The wave functions of the A-dependent nonlinear oscillator are 

X 

'^m,n{Zx, y) = Zm{Zx) Yn{y) , Zx = 



with energies given by 

1 11 
em,n = Mm + t'n = {[m + + {n + 1 - -(m + n)A 

So the total energy E„i^n = (h a)e„i.„ is a linear function of A and depends, as in the A = case, 
of the sum m + n of the two quantum numbers. It is clear that ^rn,n is well defined for any value 
of A and that the following limit is satisfied 

limA^o^!m,n{zx,y) = Hm{x)Hr^{y) e" (1/2) {^''W) ^ m,n = Q,l,2, . . . 

We recall that this approach has considered the space L2(]R2,d/XA) as the appropriate Hilbert 

space. Thus, if ^rn,n{x,y) and ^r,six,y) are to wave functions representing states of the nonlinear 
oscillators with quantum number (m, n) and (r, s) respectively, then the scalar product is given by 



(*m,n , *r-,s)A = J '^m,ni^, V) '^r,s{x, y) dflA . 
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The point is that making use of the equahty 

l + Ar^ = {l + Kzl){l + ky^) 
the measure becomes as follows 



dzx dy 



in coordinates {zx,y)- Hence, making use of the factorization of the measure d/XA in coordinates 
{zx,y), wc can factorizc the scalar product in IR^ as a product of two one-dimensional scalar 
products and arrive to the following important property 



dzx 



The following two points summarize the main characteristics of the energies of the bound 
states: 

1. Spherical A < case: 

The Hamiltonian H{A) describes a quantum oscillator on the sphere {k, > 0). The 
oscillator possesses a countable infinite set of bound states ^„,m, with n, m = 0, 1, 2, . . . 
and the energy spectrum is unbounded, not equidistant and with a difference between 
the levels that increases with N 

Co < ei < 62 < 63 < . . . < e,v < Cn+I < ■ ■ ■ 
ejv+1 - ejv = 1 + (Ar+ 1) |A| , N^m + n. 

The oscillations of the wave functions are reinforced and the values of the energies En^m 
are higher than in the Euclidean A = case; that is, En,m{^) > -£'n,m(0). 

2. Hyperbolic A > case: 

The Hamiltonian H[A) describes a quantum oscillator on the hyperbolic plane 
{k, < 0). The oscillator possesses only a finite number of bound states ^n,m, with 
n + m — 0,1,2, . . . , N\, N\ < 1/A — 1/2, and the energy spectrum is bounded, not 
equidistant and with a difference between the levels that decreases with N 

60 < ei < 62 < 63 < . . . < ejVA 

Cat+i - Bat = 1 - (A^+ 1) A, N^m + n. 

The oscillations of the wave functions are smoothed down and the values of the energies 
En,m is lower than in the Euclidean A = case; that is < En^m{0). 
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The degeneracy of the energy levels is the same that in the Euclidean case. 

Figures VI and VII show the values of the energy e^.n as a function of iV, iV = m + n, for 
several values of A. Figure VI shows as < -^^^^(O) when A > (hyperbolic A > 

case) and £'„^^(A) > £'„„j(0) when A < (spherical A < case). Figure VII shows the plot 
of the energy E^.n for three different values of A > 0; it is clear that when A decreases the 
maximum of the curve moves into the up right and the number of bound sates increases. 

6 Resolution of the Schrodinger equation II 

The second alternative way of solving the Schrodinger equation (11) is using the property 
of separability in coordinates {x,Zy) with Zy defined by Zy = y/^/l + Ax^ . This second 
approach is symmetric to the first one so the solution can be directly given as 

with 

Xn{x, m, A) = rtn^x, m, A) (1 + A x^Y (i-'"^)/(2A) ^ 
Z^{zy, A) = Hmizy, A) (1 + A 4)- V(2A) . 

The interesting property of this solution is that it satisfies 

{Hi- ^n,mix, Zy) = Vn ^n,m{x, Zy) , 

H2 ^n,m{x, Zy) = Hm ^n,m{x, Zy) . 

It is interesting to relate the existence of this second alternative approach with the above 
mentioned property of superintegrability. 

There is not presently a satisfactory definition of quantum integrability. The most direct 
way of considering this question is by direct translation of the classical notions, so that, 
according to this approach, a quantum Hamiltonian H would be integrable if there exists 
a set {Ai} of n independent observables (including the Hamiltonian itself) that pairwise 
commute. If there exists an additional set of independent observables {Br} commuting with 
H then the system would be superintegrable [60]-[71]; of course the {Br} do not necessary 
commute between them and every Br only commute with some of the {Ai}. The main 
problem is the definition of independence for quantum operators since the commutation 
relation [A^, A^] = can be considered as determining a functional dependence between Aj 
and Ak. Several criteria have been analyzed; the most simple is to consider the operators {A^} 
as independent if they are obtained by quantizating classical functions which are functionally 
independent. 

On the other hand most of classical superintegrable systems are actually super separable 
(they admit separations of variables in at least two coordinate systems) and the classical 
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integrals of motion are at most quadratic in the momenta; in this case all the operators 
{Ai, Br} are first or second order differential operators and this restriction makes easier the 
study of independence (when third order operators are considered then quantum integrability 
can lead to properties rather different to those of the classical system [62, 66]). Let us mention 
that although quantum superintegrability and exact solvability are defined in different ways, 
it has been conjectured [66] that all maximally superintegrable systems are exactly solvable. 

A consequence of the quantum superintegrability is that, as we have a total of 2n — 1 
operators commuting with H, we can construct different complete sets of n commuting 
observables and, therefore, different ways of characterizing the wave function ^. Concerning 
this quantum A-dependent oscillator, it is endowed with the following three sets of commuting 
observables 

{h;,h2-\.P}, {i^-Ap,:^}, 

that correspond to three different ways of decomposing H as a sum of two commuting 
observables and also to three coordinate systems separating the Schrodinger equation 

(zx,y), {x,Zy), {r,(l)), 
and to three alternative ways of representing the wave function. 

7 Final comments and outlook 

The harmonic oscillator is not a specific or special characteristic of the Euclidean space but it 
is well defined in all the three spaces of constant curvature. In fact if we use the curvature k, 
(or A) as a parameter then we can say that there are not three different harmonic oscillators 
but only one that is defined, at the same time, in the three manifolds. This property, that 
was known at the classical level, is also true for the quantum system. 

If we consider the spherical and hyperbolic systems as a deformation of the well known 
Euclidean system (in the sense discussed in the Introduction) then this deformation appears 
as clearly asymmetric. This is a natural result since the sphere and the hyperbolic 
plane H'^ are geometrically different and, because of this, some dynamical properties, as 
the characteristics of the wave functions ^m,n ^-nd the energies Em,n, ^-Iso show differences 
depending of the sign of A. 

We finalize with some open questions. 

Firstly, we have focussed our study on the systems {x,Zy) and {zx,y) because of their 
relations with the Cartesian coordinates and the Hermite polynomials in the Euclidean 
A = case. Nevertheless the resolution in polar coordinate must also be studied. 

Secondly, we have quantized the system by analyzing the symmetries of the metric, obtain- 
ing an invariant measure and expressing the Hamiltonian as a function of the Noether mo- 
menta. The use of this quantization procedure for other systems with a position-dependent 
mass is a matter to be studied. 
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Thirdly, one of the byproducts of this study is the existence of A-dependent deformations 
(or generahzations) of the Hermite polynomials endowed with the appropriate properties 
of orthogonality, Rodrigues formula and generating functions. The first family was already 
obtained in Ref. [1] but the second family is new; it is clear that the three-dimensional 
oscillator will lead to a new third family. These polynomials are interesting and deserve a 
more detailed study not only for the quantum problem but also from a mathematical view- 
point. At this point we also recall the existence of another family of "relativistic Hermite 
polynomials" obtained by Aldaya ct al in Ref. [74, 75, 76] in the study of the relativistic 
quantum harmonic oscillator. Wc also note that the equations (15) and (18) are of hypergeo- 
metric type and the relation of the solutions with the hypergeometric series is an interesting 
problem. 

Four, it was proved that the one- dimensional A-oscillator is solvable [1] by the use of the 
Schrodinger factorization formalism in terms of first order differential operators A and 
(these operators are known as intertwining operators) . The factorization of two-dimensional 
systems still remains as a very difficult problem but, in any case, the existence of appropriate 
operators A and for this two-dimensional oscillator must be studied. 

Finally, the technique of introducing the curvature k as a parameter for the joint analysis 
of the dynamics in the three manifolds (S*^, -E^, H^) has been generalized to the Cayley-Klcin 
geometries (see Refs. [77]-[80]). In this more general case this technique is used with two 
parameters, Ki and K2, which correspond to a space M^^^^^ with constant curvature Ki and 
signature (-1-1, ^2). This formalism is more general and includes the K-dependcnt formalism 
as the particular case Ki = k and K2 = 1. Wc think that the study presented in this article for 
the quantum harmonic oscillator on the three classical spaces of constant curvature can be 
extended to the de-Sitter, Minkowski and anti-de-Sitter spaces by using this two parameters 
formalism. 

Appendix: Geodesic polar coordinates 

A two-dimensional manifold M can be described by using different coordinate systems. If we 
consider it as an imbedded submanifold of H^, then the points of M can be characterized by 
the three external coordinates, as (x, y, z) or (r, 0, 6), plus an additional constraint relation. 
Nevertheless, in differential geometric terms, a more appropriate approach is to develop the 
study by using two-dimensional systems of coordinates adapted to M. 

On any general two-dimensional Riemannian space, not necessarily of constant curvature, 
there are two distinguished types of local coordinate systems: "geodesic parallel" and "geo- 
desic polar" coordinates. They reduce to the familiar Cartesian and polar coordinates on 
the Euclidean plane and both are based on a origin point O and an oriented geodesic g\ 
through O. 

For any point P in some suitable neighbourhood of O, there is a unique geodesic g joining 
P with O. The (geodesic) polar coordinates {R, $) of P, relative to the origin O and the 
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positive geodesic ray of gi, are the (positive) distance R between P and O measured along g, 
and the angle $ between g and the positive ray gi, measured around O. These coordinates are 
defined in a neighbourhood of O not extending beyond the cut locus of O; polar coordinates 
are singular at O, and $ is discontinuous on the positive ray of gi. 

In the case of M being a space of constant curvature k, the expression for the differential 
element of distance ds"^ is given by 

dsl = dR' + sl(R) d^' , 

so that we get ds^ = dr"^ + dcj)^ for the particular «; = Euclidean case. 
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ure Captions 



Figure I. Plot of the one-dimensional potential V(X) = (1/2) (q;^x^)/(1 + Ax^), as 
a function of x, for a = 1 and A < 0. 

Figure II. Plot of the one-dimensional potential V{X) = (1/2) {a'^x'^)/{l + Xx'^), as 
a function of x, for a = 1 and A > 0. 

Figure III. Plot of the potential U^ir), a = 1, as a function of r, for k — —1 (lower 
curve), K — Q (dash line), and k — 1 (upper curve). 

Figure IV. Plot of the A-dependent Hermite function ^2(2:0;, A) as a function of Zx 
for A = (dashed curve) and A = —0.15 and A = —0.30. For very small values of 
|A| the figure is very similar to the standard Hermite curve and when the value of |A| 
increases the oscillations become stronger. 

Figure V. Plot of the A-dependent Hermite function Z2{zx.,-S) as a function of 
for A = (dashed curve) and A = 0.15 and A = 0.30. In this case when the value of A 
increases the oscillations become softer and for A > 0.5 the Hermite function becomes 
not normalizable. 

Figure VI. Plot of the energy cat as a function oi N , N ^ m + n, ior h. ^ 0.30 

(lower curve) and A = —0.30 (upper curve). The thick points (A^, cat), corresponding 
to the values = 0, 1,2, represent the three bound states existing for A = 0.30 and 
the first three bound states for A = —0.30. The straight line (dashed line) placed in 
the middle corresponds to the linear harmonic oscillator. 

Figure VII. Plot of the energy e^.n as a function of A^, A^ = m + n, for A = 0.45 
(lower curve), A = 0.30 (middle curve) and A = 0.15 (upper curve). The curves 
also show the plot of the points (N^cn) for the values A^ = 0, 1, A^ = 0,1,2, and 
A^ = 0,l,...,6, respectively. Every thick point represents a certain number of bound 
states with the same energy i?m,„ and characterized by quantum numbers m and n 
such that n + m = N. When A decreases the maximum of the curve moves into the up 
right, the number of bound sates goes up and in the limit A — > the curve converges 
into a straight hue parallel to the diagonal (dashed hue) . 
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Figure I. Plot of the one- dimensional potential \^(A) — (1/2) (q;^x^)/(1 + A 
function of a;, for a = 1 and A < 0. 
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Figure III. Plot of the potential U^ir), a = 1, as a function of r, for k = —1 (lower 
curve), K — (dash line), and k — 1 (upper curve). 




z 



Figure IV. Plot of the A-dependent Hermite function Z2{zx,A) as a function of Zx for 
A = (dashed curve) and A = —0.15 and A = —0.30. For very small values of |A| the 
figure is very similar to the standard Hermite curve and when the value of |A| increases the 
oscillations become stronger. 
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Figure V. Plot of the A-dependent Hermite function Z2{zx,A) as a function of for 
A = (dashed curve) and A = 0.15 and A = 0.30. In this case when the value of A 
increases the oscillations become softer and for A > 0.5 the Hermite function becomes not 
normalizable. 
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Figure VI. Plot of the energy cat as a function of A^, = m + n, for A = 0.30 (lower 
curve) and A = —0.30 (upper curve). The thick points (A^, e^v), corresponding to the values 
A^ = 0, 1, 2, represent the three bound states existing for A = 0.30 and the first three bound 
states for A = —0.30. The straight line (dashed line) placed in the middle corresponds to 
the linear harmonic oscillator. 
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Figure VII. Plot of the energy Efn,n as a function of N, N = m + n, for A = 0.45 (lower 
curve), A = 0.30 (middle curve) and A = 0.15 (upper curve). The curves also show the plot 
of the points {N, En) for the values = 0, 1, = 0, 1, 2, and A" = 0, 1, . . . , 6, respectively. 
Every thick point represents a certain number of bound states with the same energy Ejn,n 
and characterized by quantum numbers m and n such that n + m = N. When A decreases 
the maximum of the curve moves into the up right, the number of bound sates goes up and 
in the limit A ^ the curve converges into a straight line parallel to the diagonal (dashed 
hne) . 
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